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ABSTRACT 


An experimental investigation has been conducted to 
examine the convective heat transfer in straight and curved 
ducts of rectangular cross section. The experimental con- 
figuration was modeled as infinite parallel plates with one 
wall at a constant heat flux and the opposite wall adia- 
batic. Experiments were conducted =: steady state for 
turbulent flow using average Nusselt numbers to compare the 
heat transfer characteristics of the straight and curved 
sections. 

The development of Taylor-Gortler vortices in the curved 
section proved to enhance the heat transfer rate over that 
of the straight section. Improved heat exchanger and turbine 
blade design could both result from a better understanding 
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A.  TAYLOR-GORTLER VORTICES 

Since the early part of the twentieth century, consider- 
able research has shown that fully developed laminar flow 
along a concave wall does not remain two-dimensional [Refs. 

1, 2, 3]. The flow instead forms a system of spiral vorte#ees, 
of counter rotating pairs, whose axes are aligned in the 
principle flow direction. This phenomenon is the result of 
the variations in the centrifugal forces acting on the 

fluid particles, and is known as Taylor-Gortler vortices. 
Figure 1 illustrates the type of fluid motion just 

described. 

In a channel that is curved in the streamwise direction, 
those fluid particles located in the center of the flow cross- 
section are subjected to higher centrifugal forces than those 
fluid particles traveling along the channel's boundary wall. 
As a result, the fluid in the center of the channel moves out- 
wardly towards the concave boundary. As the process continues, 
the fluid particles near the boundary wall, move in a spanwise 
direction, and finally radially inward replacing the outwardly 
moving particles. These particles then come under the same 
centrifugal forces and the process repeats itself. The re- 
sulting cyclic motion causes Ehe Eermationge: 76 cCoumem 
rotating Taylor-Gortler vortices, considered to b2 primar y 


a laminar flow phenomenon. 
14 
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Primary Flow Direction 






po ' 


(us 


ele jal 
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Channel. 
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It has been observed [Ref. 4] that the heat transfer meee 
from flow along a concave curved wall is greater than that 
for flow along a straight wall, a phenomenon attributed to 
the additional mixing provided by the secondary motion of 
the Taylor-Gortler vortices. It 1s thought thatthe emo oo 
hatching observed on reentry vehicles can be explained in 
part by the presence of streamwise vortices similar to the 
Taylor=Gortler yvortice pheneme ten. 

There are many possible applications that could result 
from a more thorough understanding of Taylor-Gortler vortices 
and their effect on heat transfer and fluid flow character- 
1stics. Two such applications could lead to improved heat 


exchanger designs and improved turbine blade cooling [Refs. 


Eee bo OR 

The instability of an in iseid £luid £1 owrn@. case 
curved boundary, was first cc .nsidered by Lord Rayleigh in 
1916 [Ref. 8]. By assuming that the fluid was non-viscous, 
he determined that for the motion of an inviscid Fluid te 
remain stable, its circulation must increase with increasing 
radius. G. I. Taylor, in 1923, [Refts. L, 9) 9eecontinved se aaee 
study with an extensive analytical and experimental study 
of viscous fluids. His investigations focused on the flow 
between two cylinders, in which the inner cylinder rotated 


while the outer cylinder remained stationary. Taylor 


16 


ascertained that such flows become unstable when the value 


of a dimensionless parameter exceeded a critical value of 





41.3. The parameter, known as the Taylor number is defined 
as: 
d 
pe Roy 


where ‘d' is the width of the gap, assumed small when compared 
Bemeri: , the radius of the inner cylinder, and 'Re' is the 
Reynolds number based on the peripheral velocity of the inner 
cylinder. Taylor determined that for those cases in which 

the value of the Taylor number exceeded the critical value, 

a secondary motion developed and the Taylor vortices formed. 
Figure 2 illustrates this fluid motion. 

Instability of a similar nature is also observed when a 
viscous fluid flows in a curved channel due to a pressure 
gradient acting along the channel wall. This problem was 
Preest considered analytically by W. R. Dean [Ref. 10] in 1928, 
for a channel formed by two concentric cylinders, where the 
radius of the inner cylinder was large in comparison to the 
Spacing between the inner and outer cylinder walls. Dean 
concluded that there would be an initiation of flow instabil- 
ity, and the formation of vortices, when a dimensionless 
parameter, the Dean number, exceeded a value of 36. The 


Dean number is defined as: 


a 
De = Re SD 





NY 
\\ | 
! MK NANG 


INK 


where 'd' represents the channel half-width, 'Ri' is the 
inner cylinder radius, and 'Re' is the Reynolds number based 
on the mean velocity of the undisturbed flow and the channel 
half width ‘d'. The analytical work of Dean was later veri- 
fied by W. H. Reid [Ref. 11], using an approximate solution. 

In 1940, H. Gortler [Ref. 2] studied the stability of 
laminar boundary layer profiles on curved walls under the 
influence of small disturbances. He found that these dis- 
turbances were similar to the vortices studied by G. I. 
Taylor. Using approximate numerical calculations, Gortler 
concluded that the disturbances or vortices were produced 
only on the concave boundary walls and that the overall flow 
profile appeared to remain laminar in nature. These results 
were verified, with an exact solution, by G. Hammerlin, and 
meee cad OY Hemochitem@ting [Ref.=i2])> im Y955>" A. M. Oy Smith 
completed an even more extensive numerical analysis that fur- 
ther substantiated these findings [Ref. 3]. The results of 
these numerical solutions have recently been demonstrated 
with the use of hot wire anemometry, laser doppler systems, 
and flow visualization techniques, [Refs. 13, 14]; and in 
1976, Y. Aihara [Ref. 15] conducted a non-linear analysis 
of Gortler vortices. 

As interest began to develop concerning the effects of 
the secondary flows associated with the Taylor-Gortler 
vortices, many studies were published concerning the in- 


fluences of these vortices on heat transfer in the laminar 


ro 


and turbulent flow regimes. in 1955, F. Krereh [Ret lew 
studied the influence of heat transfer with respect to the 
curvature of the boundary wall for fully turbulent tilewee 
He concluded that the heat transfer from the heated concave 
boundary wall was considerably higher than that transferred 
from the convex boundary wall of the same curvature and 
under similar turbulent conditions. 

In 1965, L. Persen [Ref. 17], considered the special cases 
of very “igh and very low Prandtl number fluids and related 
the ane se in heat transfer from a curved wall to the pre- 
sence of the Taylor-Gortler vortices. 

There has been only a limited amount of published litera- 
ture dealing with the flow and heat transfer in curved channels 
of rectangular cross-section and large aspect ratios, where 
the aspect ratio is equal to the spanwise distance of the 
channel divided by the channel height. Much of what has been 
published involves the development of numerical approximations 
and solutions for heat and mass transfer in curved ducts of 
various geometries. kK. Cheng and M. Akiyama [Ref. 18] devel- 
oped a numerical solution for forced convection heat transfer 
with laminar flows in curved channels of rectangular cross- 
section, but only for small aspect ratios. 

In 1976, A. A. Shibani and M. N. Ozisik [Ret s 19), wWemme 
matched asymptotic expansion technigues for a wide range of 
Prandtl n. ».ers, solved the heat transfer problem between paral- 
lel plates with turbulent flow, for the case of uniform wall 


temperature. 
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Other experimental and analytical studies that are worth 
citing with regard to this present study follow. Y. Mori 
[Ref. 20], obtained results for hydrodynamically fully devel- 
oped flows with constant wall heat flux, in curved channels 
Of square cross-section. W. M. Kays and E. Y. Leung [Ref. 21] 
reported solutions for turbulent flow heat transfer in a con- 
centric circular tube annulus with a fully developed velocity 
profile and constant heat rate per unit length, for a fluid 
of Prandtl number 0.7. Results for large aspect ratio 
channels with rectangular cross-sections for laminar flows 
were reported by M. Durao [Ref. 22! and J. Ballard [Ref. 23]. 
R. Holihan, Jr. [Ref. 24] reported results for laminar and 
transition flows and S. F. Daughety [Ref. 25] reported 
results for turbulent flows. P. F. Brinich and R. W. Graham 
[Ref. 26] reported results for turbulent flows in a rectan- 
gular curved channel with an aspect ratio of 6, for the 
inner wall heated, the outer wall heated, and both walls 


heated. 
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Tl. INTENT oe Sheol 


The purpose of this investigation was to examine the 
effect of streamwise curvature on the heat transfer rate in 
a curved rectangular duct of large aspect ratio, again defin- 
ing aspect ratio as the spanwise distance divided by the 
channel height. The flow velocities examined were in the 
transitior and turbulent regime. The results of this study 
were compared to the heat transfer rates for the same range 
of velocities, in a straight duct of identical aspect ratio. 

Although Taylor-Gortler vortices are considered to be a 
laminar flow phenomenon, enhanced heat transfer in curved 
channels for transition and turbulent flows has been observed. 
Even though it has not been proven that the vortices continue 
to propagate at the higher flow velocities, this enhancement 
in the transfer of heat has been attributed to the secondary 
Flow velocity components of the Taylor-Gortler vortices. It 
is believed that these secondary components transport the 
heated fluid from the outer concave wall, inward toward the 
opposite wall of the channel, displacing the cooler fluid 
particles and causing them to move toward the heated concave 
wall. It was expected that similar results would be observed 
in -thas study- 

This investigation was conducted using a single channel 


w.th a rectangular cross-sectio.1 and constant aspect ratio. 


The channel incorporated both a straight test section and a 
oe 


curved test section. The results obtained in each test section 
at approximately the same flow rates were compared in an effort 
to determine the effects of the Taylor-Gortler vortices on the 
transfer of heat. Also, the results of both the straight and 
curved sections were compared to the results of Ballard [Ref. 
23), Holihan [{Ref. 24], and Daughety [Ref. 25]. 

The straight section results of this study were compared 
to the Dittus-Boelter equaticn for flows in circular tubes 
[Ref. 26], the results of Kays and Leung for turbulent and 
transition flow in annular passages [Ref. 21], and the 
analytical solutions cf Petukhov and Popov for transition 
flows in circular tubes [Ref. 27]. 

The curved section results of this study were compared 
with the results of Brinich and Graham [Ref. 28] for turbulent 
flow in a rectangular curved channel, and also to the solutions 


of Petukhov and Popov, mentioned above. 
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IIIT. EXPERIMENTAL WORK 


A. DESCRIPTION OF THE APPARATUS 

To meet the objectives of this investigation, a channel 
of rectangular cross-section was used. Details of its con- 
struction are described in references 22 and 23. Plexiglas, 
0.635 centimeters thick, formed the inner and outer Wale 
the channel. These walls were separated by 0.635 centimeter 
Spacers that also served as the sides of the channel. 

The channel, shown in Figure 3, had a straight section, 
122.0 centimeters in length, followed by a curved section of 
180 degrees of arc, followed by another straight section, 91.4 
centimeters in length. The inner curved wall of the channel 
had a radius of curvature of 30.5 centimeters. The rectan- 


gular channel was 0.635 centimeters high and 25.4 centimeters 


wide resulting in an aspect ratio of 40. The hydraulic dia- 
meter was 1.229 centimeters. The cross-sectional area of 

the channel was 16.13 square centimeters. The wetted peri- 
meter was 52.07 centimeters. Figure 4 shows a cross-sectional 


view representing both the straight and curved sections of 
the channel. 

An entrance bell constructed of plexiglas was connected 
to the straight section of the channel. It was designed and 
manufactured according to ASME nozzle standards, with an 


elliptical curved base on a major axis equal to ten inches 


24 
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and a minor axis of one inch. Cheese-cloth was attached to 
the entrance nozzle to prevent foreign matter from entering 
the flow channel. An aluminum exhaust nozzle was attached 

to the exit of the channel and directed the flow from the 
channel into two inch diameter pvc piping. The two inch pve 
piping contained a measuring section which could be fitted 
with standard ASME concentric orifices. In this study two 
orifice plates were used. An orifice with a diameter of 
0.5325 inches was used for hydraulic Reynold's numbers 

below 5000, and an orifice with a diameter of 1.065 inches 

was used for hydraulic Reynold's numbers above 5000. Pressure 
taps on either side of the orifice, at 1 and 1/2 pipe diameters, 
and connected to manometers were used to calculate the mass 
flow rate of the working fluid, which was air at room 
temperature. 

The flow of air was drawn through the channel by an 
electrically driven Spencer Turbo Compressor, rated 30 horse- 
power at 3500 rpm, and 550 cubic feet per minute at 70 degrees 
fahrenheit and one atmosphere. 

The channel contained two test sections from which experi- 
mental heat transfer data could be obtained. The straight 
test section was located 120 centimeters downstream of the 
entrance bell to ensure hydrodynamically fully developed 
flow. The straight test section was 29.2 centimeters in 
length with a heated test section area of 741.7 square centi- 


meters. The curved test section was located in the lower 
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half of the curved portion of the channel. It was 28.3 
centimeters in length, subtending an arc of 53.1 degrees, 
and having a heated test section area of 718.4 square 
centimeters. 

In each of the heated test sections, Temsheet, a carbon 
impregnated porous paper with a uniform electrical resist- 
ivity, was glued to the outer wall. Joulean heating was 
used to heat the flow of air through the channel. Since the 
electrical resistance of the Temsheet is not constant, but 
varies slowly with temperature, a precision resistor with an 
electrical resistance of 2.0262 ohms was connected in series 
with the Temsheet to allow the calculation of the instantan- 
eous power being supplied. Detailed schematics of each test 
section are shown in Figures 5 and 6. 

The variables that were measured and used in this invest- 
igation were: 

(1) the temperature of the air entering the channel (T._) 


in 
) 


(2) the temperature of the air leaving the channel (ie 
(3) the temperature of the heated outer wall for each 
test section (T_) 
WO 
(4) the temperature of the unheated inner wall for each 
test section (T_.) 
Wi 
(5) the temperature between the three layers of insulation 
at each test section (T.__) 
ins 


(6) the voltage across the precision resistor (V_) 


(7) the voltage across the Temsheet or heater (Vi) 
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Figure 6. Curved Test Section, Detailed Schematic. 
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(8) the temperature of the air at the orifice OTe) 

(9) the pressure upstream of the orifice (Pup) 

(10) the difference in the upstream and the downstream 
pressure across the orifice (AP) 


) 


(ll) the atmospheric pressure areca 

A total of sixty-one Copper-Constantan, 30 gauge thermo- 
couples were located throughout the test apparatus to record 
the desired temperatures. Readings from four thermocouples 
Spanning the width of the channel were averaged to give the 
temperature at the entrance and exit to the channel. Five 
thermocouples were connected in parallel and inserted between 
each of the layers of insulation at each of the test sections. 
In addition, individual thermocouples were positioned in each 
of the test sections, eight in direct contact with the Temsheet 
and eight in the plexiglas of the unheated inner wall. The 
beads of the thermocouples in contact with the Temsheet were 
electrically insulated with ENMAR Heat Resisting Glyceryl 
Phthalate. Small diameter holes were Crilled in the plexi- 
glas to allow the installation of the thermocouples. One 
thermocouple was inserted in the two inch pvc pipe cownstream 
of the orifice in accordance with the ASME Power Test Code 
to record the temperature of the air flowing through the 
Oritice, 

The thermal insulation consisted of 1/2 inch layers of 
Armstrong Armaflex 22 Sheet Insulation, a flexible foamed 


plastic material. By positioning thermocouples between 


ol 


these insulation Il¥yers, the eat lost to the environment 
could be computed t each tes. section. The heat loss 
through insulation surrounding the rest of the channel due 
to the heated air flowing through the channel was small in 
comparison. See Figure 7 for a detailed sketch of the 
thermocouple placement and attachment. The insulation was 
held in place by the use of velcro strips attached to the 
frame. The entire channel and all connections between the 
channel, the pve piping, and the orifice were sealed with 
General Electrical Silicone Rubber Sealant Caulk to ensure 
there was no leakage of air into the channel or the piping 
which could affect the recorded temperatures or mass flow 
eae Cs 

The data acquisition system used for this experiment 
was a Hewlett Packard 3054A Automatic Data Acquisition/ 
Control System consisting of a 3456A Digital Voltmeter and 
a 3497A Data Acquisitien/Control Unit. Also used 1njcan. 
Junction with this data acquisition system, were a Hewlett 
Packard 9326 Computer terminal and 2671G Printer. Pressure 
measurements were taken from Meriam vertical manometers, one 
whose fluid was water with a 0 to 60 inch range, another 
using water with a 0 to 30 inch range, and a third whose 
fluid was mercury, calibrated to read inches of water with 
a range of 0 to 415 inches. A photceqrapliet eneoscrannel 
and associated st equipment are shown in Figures 8, 9, 


and 10. 
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Photograph of Test Channel. 
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Photograph of Flow Measuring Apparatus. 
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Sheets of cholesteric liguid crystal slurries with a 
temperature range of 40-43°C were used to provide gualita- 
tive measurements of the surface temperature distribution 
on the Temsheet. The liguid crystals exhibit dramatic color 
changes over a small temperature differential and allow a 
reasonable comparison to be made between the heat transfer 
processes on each test section by observation of the color 
pattern displayed for varying temperature gradients and flow 
velocities. Specifically, the occurrence of multi-colored 
Stripes aligned in the flow direction on the curved wall 
would indicate regions of high and low temperatures and con- 
Firm the presence of Taylor-Gortler vortices. 

The addition of this liguid crystal sheet necessitated 
Emealeulation Of its Overall heat transfer coefficient. To 
determine the true surface temperature of the heated wall 
the temperature drop across the sheet must be known. The 
sheet was composed of 4.5 mils of a liguid crystal slurry 
with approximately the same thermal conductivity as body 
Rios Of acrylic backing, and a 2 mil mylar cover; 
resulting in an overall heat transfer coefficient of 
522 w/°C. Insulation opposite the heated sections was 
fitted with removable panels to allow viewing and photo- 


Seapning of the liquid crystals. 


B. EXPERIMENTAL PROCEDURES 
The experimental procedures followed were the same for 
the straight and curved test sections. The results obtained 
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from the straight test section served as a baseline for the 
comparison of the curved section results. Experiments were 
conducted by varying the mass flow rate of air corresponding 
to Reynolds numbers from approximately 800 to 25000, where 
the Reynolds numbers are based on hydraulic diameter. For 
each mass flow rate, data was taken automatically by the data 
acguisition system and the data was immediately reduced and 
printed. Preliminary runs were performed to determine the 
time reguired for the test rig to come to steady state. 
For this experiment three hours were allowed for steady 
state to be reached. After three hours, data was taken at 
ten minute intervals to ensure that a steady state condition 
had been achieved. The criteria for steady state was based 
on three variables: 

(1) the difference between inlet and outlet temperature 

(2) the heated boundary wall temperature 

(3) the unaccounted neat loss in the system 
When the first two of these variables varied by less than 
two percent over a ten minute interval and the third began 
to fluctuate, it was considered that a steady state condition 
had been achieved. 

It was also necessary to determine the approximate power 
setting required to bring the heated wall of each test section 
tc approximately 50 degrees Centigrade. It had been deter- 
mined by Holihan [Ref. 24] that this temperature was suffi- 
cient to ensure a 20 degree Centigrade difference in the 


heated wall and the unheated wall. 
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For each run, the atmospheric pressure, the pressure 
difference across the orifice, and the pressure upstream of 
the orifice were entered into the computer for the mass flow 
rate calculations. All other data was acquired by the data 
acquisition system directly from the experimental apparatus. 

To compute the instantaneous power supplied the following 


relationship was used: 


The precision resistor voltage and the heater voltage were 
read by the data acquisition system and the resistance of the 
precision resistor was a known constant. A diagram of this 
@meult 1S shown in Figure ll. 

Upon completion of selected data runs photographs were 
taken of the liquid crystals on the heated section to provide 
emevisual record of the airflow patterns. 

For Reynolds numbers ranging from approximately 2400 to 
5600, it was necessary to make two sets of data runs for the 
straight test section. The first, with conditions the same 
as for other flow rates, which exhibited a pronounced devia- 
tion from the near linear results; and a second, where the 
flow just prior to the straight test section was tripped 
with sandpaper, which reestablished the near linear results. 
This reproducible anamoly indicated the presence of a tran- 
Sition region between laminar and turbulent flow. 
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Figure ll. Circuit Used to Measure Power Supplied to 
the Temsheet. 
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IV. PRESENTATION OF DATA 


A. ANALYSIS 

The constant heat flux surface for the heated wall of each 
test section was provided by the uniform electrical resis- 
tivity of the Temsheet. The insulated unheated wall was con- 
Sidered adiabatic since the heat losses through that wall were 
negligible. The channel's large aspect ratio of 40 allowed 
the modeling of the channel as parallel plates. The straight 
portion of the channel upstream of the straight test section 
was of sufficient length to ensure that the flow was hydro- 
dynamically developed for the flow velocities of this study. 
The straight portion of the channel downstream of the curved 
test section was also of sufficient length to ensure that the 
flow exiting the channel was thoroughly mixed and that the 


mmectature, T was an average bulk temperature. 


Outs 
Based on these assumptions, the experimental configura- 
tion was modeled as forced convection between parallel plates, 
for hydrodynamically developed and thermally developing flows, 
subjected to a constant heat flux. The boundary conditions 

are one wall at a constant heat flux and the other wall 
insulated. 

To analyze this problem, several quantities were defined 
as follows; 

The heat convected to the air was calculated using the 


eguation: 
41 


QO. =mc in = eee) 


al Ws DABYr. Olte an 
where Cpair was the specific heat of the air at constant 
pressure and 'm' was the mass flow rate of the air. Mass 


flow rate was calculated using standard ASME flow measure- 
ment technigues [Ref. 33]. 
The average heat transfer coefficient between the heated 


wall and the flow of air in the channel was def .ed by the 


equation: Q 
ae ghia 
Apr, a 
where Gi 1s defined above, ‘Apr, was the area of the 


Temsheet in the test section, and 'AT' was the difference 


between the arithmetic average heated wall temperature (i 


) 


O 


and the average bulk air temperature (Thay): The average 
bulk temperature was defined as the arithmetic mean of the 


he 


entrance and exit temperatures (T._, T 
in Oue 


The average Nusselt number was then calculated using: 
aD 
Ve, = ue 


hd Kae 
ellis 


ieee ae Sea Walieaom ‘Dug’ is the hydraulic diameter and ‘koe 
is the thermal conductivity of air. 
The Reynolds number was calculated for each test run as 


follows: 


hd 


CG "aie 


| 
= 


hd 
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where again 'm' and 'D,.' are the mass flow rate and hydrau- 


hd 
Hic diameter of the channel, AR! 1s the cross-sectional 
area of the channel, and ee is the dynamic viscosity of 
the air. 


For the curved section runs the Dean number was defined as: 


hd 
hd Ri 


where ‘Rey is the Reynolds number based on the hydraulic 
diameter, ‘Dig’ is the hydraulic diameter, and 'Ri' is the 
radius of curvature of the unheated inner convex wall. 

A sketch of the control volume and a set of sample calcu- 
lations for one test run of the curved section are given in 
Appendix B. A sample calculation for the uncertainty analysis 


1S given in Appendix A. 


Pee RESULTS 

The data obtained from each experimental run was evaluated 
utilizing the expressions described in Part A above. The 
major parameters resulting from this evaluation are shown in 
Tables I through V. Tables I through III contain the straight 
section results and Tables IV and V contain the results of 
the curved section. A plot of the average Nusselt number 
verses Reynolds number is given in Figure 12 for the compari- 
son of the results. Figure 12 shows the flow passing through 


the transition regime of the straight section in both the 
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Table I 


Summary of Straight Test Section Results 
(Large Orifice) 


| A 
“Sha io (W7/mC) a “ha 
4973 47.09 30.80 20.61 14998 
5035 67.70 3 SS 2560 L7 6Z2e 
6326 a5 oe Seles) Zoe oo 18.2 
7096 75.40 44.41 22.88 21.54 
7989 98.82 AO Zola 23.08 
8935 82.48 52.05 21.36 25.10 
9969 101.85 56.20 DA 27 ee 
11520 ae 63.12 25.90 30.52 
ean 120.00 67.01 24.13 32.35 
14357 eo 73.52 22.40 35.43 
16137 120409 77.78 Dea6 37-50 
17990 124.29 83.17 20.14 40.55 
20629 25.3 0 69553 Sic 6 43.62 
24597 127.56 lianas 17.00 48.63 


Reynolds number and Nusslet number are based on hydraulic 


G@taemetrer-. 


Table II 


SUMMeatayeet Straight Test Section Results 
(Small Orifice - Untripped) 


“Shd sie ae es nnd 
1524 aay 14.74 24.84 7226 
2162 Beg 18.92 228 9.20 
2686 34.81 19.03 24.65 9.23 
7a 39.19 Dea 2G 10.57 
4063 AGeeei Dana re yl 11.00 
4463 39.28 Dees 2208 11.60 
4851 52.09 26.41 26.58 12.77 
5223 51.15 30.97 22.25 15.18 


Reynolds number and Nusselt number are based on hydraulic 


diameter. 
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Table III 


Summary of Straight Teste ,Seculem Results 
(Small Orifice ~- Tripped) 


Reng Cair ‘a vi it 
(W) 
TI TG 41.74 ok Sale © 8.362 
2354 34.55 lL oeor: Dood) Oo 
27s 40.11 22S Zee i Oras 
225 ao! 2A oe L944 ll 33 
3626 DO m6 7 Zoey 24423 13.72 
4639 50.47 33.464 20 FAG 16.44 
Ser Soe 39735 Len Oras 
7014 1255 45 Ao Zee 22.29) 
Reynolds number and Nusselt number are based on draulae 
diameter. 
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Table IV 


(Large Orifice) 


Cair 2 
(W) (W/mC) 
79.49 41.71 
Tos 0 44.95 
82.47 pile 5 
81.40 54.48 
OF 75 59.04 
IQ eo 64.62 
Oe 19 68.81 
MOS2 52 2a) 4 
110.84 F220 
(eld SZ TIES 
PAG en 4 Seas G 
Peo oo SGnns 9 
eS iO 
EZ 423635 Hl SO4 
AG ol 2624 3 
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Reynolds number and Nusselt number are based on 


diameter. 
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Summary of Curved Test Section Results 
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tripped and untripped condition. All further comparisons 
will be conducted using the data from the tripped flow. 

The results indicate an increase in the rate of heat 
transfer with increasing Reynolds number for both the 
Straight and curved test sections. In addition, the heat 
transfer rate was higher in the curved section for each 
Reynolds number investigated. 

A straight line least squares correlation of the data 


for the present study resulted in: 


—_ O27 
Nu = 0.02 Rev g 2300- Rey a <A EG) 


for transition flow in the straight test section; 


— 2 
Nu = 0.22 Re 0.33 


end 1500 < Re 


ae 53010 


FOr transition flow im the eurved test scecelone 


Nu = 0.07 Re 0.65 


ae 8900 < Re 


eae Zo 0010 


for turbulent flow in the straight test section; and 


O55 


Nile = Oeedles Revd 


8400 < Re < 25000 


hd 


for turbulent flow in the curved test section. 

In the short laminar regime of this study the rate of 
heat transfer was about forty-five percent higher in the 
curved section, indicating the presence of Taylor-Gortler 
vortices. As flow progressed through the transition regime 


a loss of enhancement of heat transfer rate, from forty to 


a0 





PeeluUmeWwoley percent, was observed as Elow over the straight 
section went from laminar to fully turbulent. In the low 
end of the turbulent regime, Reynolds number between 5400 
and 12000, the rate of heat transfer was about thirteen 
percent higher in the curved section. As the Reynolds 
number continued to increase beyond 12000 so did the enhance- 
ment of heat transfer, to approximately twenty-Six percent 
at 24000. This compares favorably with earlier studies by 
Kreith [{Ref. 16!, where he reported an increase in the rate 
of heat transfer along a concave wall of from twenty-five 

to sixty percent for Reynolds numbers, based on hydraulic 
diameter, between 107 and ier Holihan reported an increase 
of fifteen percent for laminar flow and thirty percent for 
transition flow [Ref. 24]. Ballard reported an increase of 
eleven percent in the heat transfer rate for laminar flows 
[Ref. 23]. Daughety reported an increase of twenty percent 
in the heat transfer for turbulent flows [Ref. 25]. 

The actual presence of Taylor-Gortler vortices was 
verified by the observation of the liquid crystal isotherm 
distribution in the curved section. At low values of 
Reynolds number and through most of the transition regime, 
liquid crystal bands oriented in the streamwise direction 
were visible (Fig. 13). As the Reynolds number increased 
turbulent forces began to dominate in the curved section, 
and less evidence could be detected for the existence of 


the vortices (Fig. 14). Over the entire range of Reynolds 


SEE 


numbers studied for the straight section, the liquid crystals 
formed essentially straight bands perpendicular to the flow 
direction, indicating approximately two-dimensional flow 
(Fig. 15). Disruption of flow in the transaerean regitom 

due to the insertion of the sandpaper is illustrated in 


Figures 16 and 17. 


2 








Figure 13. Liquid Crystal Representation of Airflow. 
(Curved Section, Reig = E00) 
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Figure 14. Liquid Crystal Representation of Airflow. 
| (Curved Section, Reng = 9400) 
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Figure] .5- 


Liquid Crystal Representation of Airflow. 


(Straight Section, Re 


54 


hd 





12700) 


Figure 16. Liquid Crystal Representation of Airflow. 


(Straight Section, Rey 4 = 4000, Untripped) 


Figure 17. Liquid Crystal Representation of Airflow. 


(Straight Section, Reg = 4000, Tripped) 
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V. DISCUSSION ANDRE CONCEUSEC iS 


Holihan [Ref. 24], in his study, determined that thease 
was a negligible difference in the fluid bulk temperature 
and the unheated wall temperature. He additionally demon- 
strated that the radiated heat transfer from the unheated 
wall surface was minimal. These results were verified by 
Daughety [Ref. 25] for the present study and provide the 
basis for the assumption that the heat transfer to the air 
Flowing through the channel was solely by convection from 
the heated Temsheet. 

As was mentioned earlier, the high aspect ratio of the 
channel provided the basis for the assumption that the ex- 
perimental apparatus, as configured, could be modeled as 
infinite parallel plates. Holihan's data [Ref. 24] also 
tended to verify this assumption. His experimental data, 
in the laminar flow region, approached the theoretical limit 
for average Nusselt number for parallel plates with one 
wall heated at a constant heat flux and the opposite wall 
adiabatic [Ref. 29]. 

Based on the assumptions mentioned above, comparisons 
were made with other analytical and experimental results 
that were of the same or similar problem. Again the problem 
being, flow between infinite parallel plates with one wall 


at a constant heat flux and the opposite wall adiabatic. 


545) 





A comparison of the experimental results of this study 
and the experimental results of Ballard [Ref. 23], Holihan 
[Ref. 24] and Daughety [Ref. 25] are shown in Figure 18, 
for the straight section, and Figure 19 for the curved 
section. Figure 20 shows a compilation of all the data on 
one plot. Ballard's data was strictly in the laminar 
region; Holihan's data covered the laminar as well as tran- 
Sition regions, and Daughety's data was strictly in the 
turbulent region. The channel aspect ratio was 40 for each 
of these studies and the exverimental procedures were similar 
in each case. Ballard and Holihan compared their data to 
the analytical studies by McCuen, [Ref. 30], for heat 
transfer between infinite parallel plates with constant 
wall temperatures and Neat flux. Daughety compared his 
Gata with the results of Shibani and Ozisik [Ref. 19], and 
Kays and Leung [Ref. 21]. Ballard's and Holihan's data 
plotted above analytical solutions while Daughety's data 
plotted below. Differences were attributed to the variations 
in the geometries used in the studies, side wall effects, 
inability to account for all of the heat transfer processes, 
and the limitation inherent in any experimental work. 

A comparison of analytical and experimental results 
with those of the present study for the straight test sec- 
tion 1s shown graphically in Figure 21. The results of 
Petuhkov and Popov [Ref. 27] are shown for turbulent heat 


meancrers IM round cubes. The correlation used was 
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N a (f/8) Re Pr 
aa (273) 
Ky (£) +k, (Pr) v£/8 (Pr ~1) 


where ‘f* the friction factor, has been calculated using the 


Filonenko equation given in Karlekar and Desmond [Ref. 32] as: 


ibe (1.782 Wieg@Re = 1.64) 7° 
and 
kK, (f) =1+41.4¢ 
Lee 
k,(Pr) = 11.7 + ——- 
2 py (1/3) 


Tne Dittus-Boelter equation [{Ref. 26] for heat transfer in 


a straight tube for constant wall temperature is plotted 


Peon: 


Nu = 0.02 Rego 


POr a Prandtl number equal to 0.71. Additionally, the 
results of Kays and Leung [Ref. 21] for heat transfer in 
annular passages are shown for r* = 1.0, which equates to 
parallel plates. When extrapolated down these results 
continue to correlate well with present data. 

The curved section data is plotted in Figure 22, and 
is compared to the experimental work of Brinich and Graham 
[Ref. 28]. Also plotted as a reference are the analytical 
Solutions Of Petukhoy and Popov [Ref. 2/7]. The accuracy Of 


the data points used from Brinich and Graham are subject to 
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REPRODUCED AT GOVERNMENT EXPENSE 
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errors, in that the actual values are not given in their 
study and had to be taken from a plot of Stanton number 
versus Reynolds number. 

Disregarding the small differences between the results 
and correlations, the trend appears consistent. Heat trans- 
fer rates can be seen to increase both with increasing 
Reynolds number and curvature of the channel. The presence 
of Taylor-Gortler vortices in the laminar and transition 
regime seems to be the major contributing factor; whereas 
increased turbulence, possibly caused by the formation of 


these vortices, seems to dominate at higher Reynolds numbers. 
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Wie RECOMMENDATIONS 


The amount of experimental work remaining in the study 
Seemeatstransfer in curved channels is significant. MThe 
present experimental apparatus has been used over a range 
of Reynolds numbers throughout the laminar, transition and 
into the turbulent regime. Additional experiments probing 
further into the turbulent regime should be conducted to 
better understand the heat transfer processes occurring 
there. This would require a more powerful compressor to 
generate the higher Reynolds numbers and a larger power 
supply to sufficiently drive the heated surfaces. Finally, 
increasing the number of thermocouples at least two-fold 
and the introduction of pressure transducers which could 
be connected directly to the data acquisition system would 
both prove to increase the accuracy and reliability of the 
experimental apparatus. 

To gain further insight into the effect of curvature 
Simeene, rate of heat transfer, additional boundary conditions 
for the test should also be considered. An apparatus possess- 
ing a heated wall on each side of the test section might ke 
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APPENDIX As EXPERIMENTAL UNCERTAINAY 


The uncertainties for the major variable in the experi- 
ments, were calculated in accordance with the method des- 
cribed by S. Kline and F. McClintoch [Ref. 29]. The estamaces 
of the uncertainty in the measured quantities were made con- 
servatively. As a result, there is considerable confidence 
in the uncertainties as calculated. 

The following eguations were used to calculate the 


uncertainties: 
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Licm— alueS obialnmea for themlumcertagmties for tie values 


in the curved section test run at a Reynolds number of 


13000 are: 

Quantity Uncertainty 
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The major source of uncertainty in the average Nusselt number 
is the uncertainty in the vressure readings from the mano- 
meter, which affects the mass flow rate calculation. [In 
order to obtain a more accurate pressure measurement, it 

has been recommended that pressure transducers be used 


instead of vertical manometers in future studies. 
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REPRODUCED AT GOVERNMENT EXPENSE 


APPENDIX B: SAMPLE CALCULATIONS 


Figure 23 pnelow, shows the major heat transfer components 
for each of the test sections. The sample calculations that 
follow, demonstrate the methods used by the computer to calcu- 
late these components as well as the Reynold's number, average 
heat transfer coefficient, and average Nusselt number for each 
set of data. The sample calculations are for the curved sec- 
tion, but those for the straight section are similar. Reynolds 
number, Nusselt number, and Dean number are based on hydraulic 


diameter. 
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SAMPLE CALCULATION DATA 


Vy = (555 28..V 
Ve = "4 ony, 
Ror = 2.0262 98 
ae = SO0e ieee. Hg 
AP =e? 7 0 CO meena HO 
Py =) 933507 in. HO 
oT = 227010 C 
6m 
Toe = 26272 C 
7 =e a 
wO 
Meee ee eae 
wi 
me = 39.67 € 
insl 
ne oo ome 
ins2 
eis = 26.57 C 
K_. = 36708 lom vai 
air 
ca moe J/Kg K 
ne -6 
Hee = (e250) x 210 Kg/ms 
K. = 3.6495 x 10°° W/mK 
ins 
Oe = 0.01270 m 
ins 
C0 = 0.70 
E. = 0.40 
wi 
9 = 5.669 x 100° wim-n 
8 = pS 
Y = 1.402 
_ Z 
ge = 1 Kg m/NS 
R = 286.8 Nm/Kg K 
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soewa. 7° 
_ Z 
Ant, = O730° m 
_ 2 
Aa ~ 0016 m 
2 
a ee = OU ZO eel 
D =i een OG Saeed 
Ee 
D8 = -0508 m 
pipe 
D ; = 0 20 om 
@isdei 
Bae = . 2 COM 


TEMPERATURE CALCULATIONS 


1. Bulk Temperature (Thik) 


Bethe) sone) 23 ol ee 6. 78. 
ie =e = sO 224.90 C 
blake 
2 Z 
2. Mean Temperature Difference (AT) 
AT = Pere Drea e = 44.74 - 24.90 = 19.85 C 
3. Temperature Difference in the Insulation (AT...) 
Te ly. ae Oey esse = 16-2956 
ins cers ins2 


POWER CALCULATIONS 


1. Power Supplied (Q,,.) 
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a. Radiation Resistance (Rp) 
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De MASS FLOW RATE CALCULATIONS 
1. Pressure Conversions 


a. P_, =30.11 in. Hg x 3374.1 = 101594.15 N/m? 





b. AP = 12.0 in. H,O x 248.64 = 2983.68 ac 


Cee.) — loomed Ooze se G4) 1OLSoteIis = 103707.59 
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Density of Air foe! 
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Expansion Factor (Y) 
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Mass Flow Rate (m) 
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Meets t Line 


Assume a Reynolds number _. = 40000 
pipe 


Obtain a value for K, the flow coefficient, from 


reference 30. R= 63350 

Solve for m. m= .0302 Kg/s 
mo Di. = 
= pipe _ (mest Z acess) 

SOEVE FOS Re pon 5 eNClO0e nmencnac ane 

pipe’ air 
Re_. = 40713 
pipe 


Check convergence and repeat process if necessary. 
(Convergence if difference less than .001) 


m= .0302 Kg/s 


REYNOLDS NUMBER CALCULATIONS 
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Pitino HEAT RANGE ER, CORPEIGCIENT CALCULATION 
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APPENDIX C: (CORRELATIONS 


The correlations obtained for the present study were 
obtained using the method of least squares for a first degree 
polynomial, as outlined by C. F. Gerald [Re&. 31]. iIngeaaee 
method the values of Reynolds number (x) and Nusselt number 


(Y) were first converted to their natural logarithmic value. 
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Next the summations of Xi, Xi, Yay and xi Y; were calculated. 


a: 


The values of these quantities were then placed in a matrix: 


and the resulting simultaneous equations 
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oy 
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The equation of the line then became: 


Ln Nu = as “++ a,in Re 


Or 


The standard error for the present study correlation was 


calculated by the following equation: 
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where: Y; 1s the actual value of the average Nusselt number 
obtained experimentally. 
Ys is the value of the average Nusselt number calculated 
from the correlation equation. 
N is the number of data points used in the correlation. 


nis the degree of the polynomial used. 
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